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ABSTRACT 

Rap1 is an essential DNA-binding factor from the 
yeast Saccharomyces cerevisiae involved in tran- 
scription and telomere maintenance. Its binding to 
DNA targets Rap1 at particular loci, and may 
optimize its ability to form functional macromolecu- 
lar assemblies. It is a modular protein, rich in large 
potentially unfolded regions, and comprising BRCT, 
Myb and RCT well-structured domains. Here, we 
present the architectures of Rap1 and a Rap1/DNA 
complex, built through a step-by-step integration of 
small angle X-ray scattering, X-ray crystallography 
and nuclear magnetic resonance data. Our results 
reveal Rap1 structural adjustment upon DNA 
binding that involves a specific orientation of the 
C-terminal (RCT) domain with regard to the DNA 
binding domain (DBD). Crystal structure of DBD in 
complex with a long DNA identifies an essential 
wrapping loop, which constrains the orientation of 
the RCT and affects Rap1 affinity to DNA. Based on 
our structural information, we propose a model for 
Rap1 assembly at telomere. 



INTRODUCTION 

Yeast Rapl (Repressor Activator Protein 1) is an 
abundant essential DNA-binding protein that plays 
multiple functional roles in vivo, and is found at pro- 
moters, silencers and telomeres (1,2). It is a key transcrip- 
tional activator of many coregulated genes (3), a 
modulator of chromatin structure at numerous yeast pro- 
moters (4), a repressor that silences transcription at HML, 
HMR and telomeres (5,6), and an essential telomeres com- 
ponent (7,8). Its association to DNA transcription sites 
affects nucleosome occupancy, and chromatin structure 
and dynamics (4,9). Rapl activates transcription at 
binding sites that can be distant by > 300 bp from the 
activated genes, in collaboration with other 
DNA-associated proteins (3,4,10). Rapl tightly and inde- 
pendently binds the double-stranded telomeric DNA with 
an average frequency of one protein every 18 bp (1 1-13). It 
recruits functional partners essential and specific for either 
negative regulation of telomere elongation, transcriptional 
repression or inhibition of non-homologous end joining 
(NHEJ) (14, 15). 

The initiation of Rapl functions relies on its inter- 
action with DNA. DNA consensus sequences have 
been published which consist of 12-14 bp (10,16,17). 
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The prototypical consensus sequence ACACCC 
RYACAYM includes two tandem half-sites ACACC or 
ACAYC and ACATY or ACAYM at bases 1-6 and 9-13 
(12,17), although thermodynamic mapping of Rapl 
binding includes three additional bases located immediate- 
ly at the 5' of this sequence (18). Rapl not only binds this 
sequence with an extreme stability (19,20), but also toler- 
ates large variation in the sequence, particularly in the 
second half-site (2,10,20). In addition to specific DNA 
binding, Rapl is able to promote modifications of the 
conformation of its target DNA site, including bending, 
untwisting and quadruplex formation (G4 DNA) 
(11,21,22). By its ability to promote single-strand 
invasion (23), Rapl appears functionally related to the 
mammalian TRF2, which is known to promote t-loop 
in vitro (24). 

Mapping studies on Rapl have identified a N- terminal 
region that includes a BRCT domain, a double-Myb DNA 
binding domain (DBD), and a regulatory C-terminal 
domain (RCT), which directly interacts with several func- 
tional partners (1) (Figure 1A). In addition to these three 
domains, 40% of Rapl peptidic chain corresponds to pre- 
dicted unstructured regions. The first 279 amino acids of 
the protein, which include the BRCT domain, can be 
deleted without affecting any known function (3), 
although yeast two-hybrid experiments have shown that 
this region is required for the interaction with the tran- 
scription factor Gcrl (25). The N-terminal and C-terminal 
parts are dispensable for chromatin opening (26), or inter- 
action with nucleosomal binding sites (9). The unstruc- 
tured linker between DBD and RCT includes a region 
required for transactivation (residues 630-695) (27) and 
a toxicity region (residues 598-616) (28). The RCT is 
required for the negative feedback loop that represses 
telomere elongation by telomerase (29,30), and for the 
establishment of a silent chromatin near telomeres 
(6,31). It interacts with the proteins Rifl, Rif2, Sir3 and 
Sir4. Rifl and Rif2 are required for the inhibition of 
telomere elongation through different mechanisms 
(15,32-35). Sir3 and Sir4 establish transcriptional 
silencing, which also requires Sir2, a conserved histone 
deacetylase (36). The binding of the Rif and Sir proteins 
to RCT appears to be mutually exclusive (37,38). Rif2 and 
Sir4 are required for NHEJ inhibition, although protec- 
tion against NHEJ is also observed with a DBD-RCT 
construct of Rapl in the absence of Rif2 and Sir4 (15). 

The presence in Rapl of several domains involved in 
specific interactions and linked by flexible regions 
suggests a high structural plasticity of this molecule as 
part of its functional competence. Recent electron micros- 
copy studies have provided the first information about 
Rapl architecture, which adopts a pseudo-ring conform- 
ation in the absence of DNA (13,39). The X-ray structure 
of the DBD in interaction with a DNA fragment is avail- 
able (12), as well as that of the C-terminal domain (38,40), 
and the NMR structure of the BRCT domain (41). 
Despite these studies, the structural determinants 
underlying the wide range of functions that Rapl is able 
to fulfill are still poorly understood. 

A crucial step to understand the role of Rapl plasticity 
is to access to its whole architecture in complex with 



DNA. The size of Rapl together with its high content of 
unstructured regions requires complementary approaches 
to characterize its 3D structure (42). Our approach inte- 
grates structural information obtained from SAXS, X-ray 
crystallography and NMR. We combined SAXS with ana- 
lytical ultracentrifugation (AUC) using three different 
constructs of Rapl: DBD-RCT, BRCT-DBD-RCT and 
full-length protein. This enabled to step-by-step build the 
architecture of the whole molecule and of its complex with 
DNA, and revealed conformational adjustment upon 
DNA binding. The crystal structure of the DBD in 
complex with a 30 bp DNA provided the structural deter- 
minants of Rapl conformational adjustment upon DNA 
binding. Finally, NMR titration of the RCT with 
full-length Rif2 and a Sir3 peptide highlighted partially 
overlapping surfaces, which remain accessible in the 
Rapl /DNA complex. 

MATERIALS AND METHODS 

Sample preparation 

The 6His-tagged Rapl[i_ 22 4], RapW-224], Rapl [117-352], 
Rapl [i_827] WT or mutant Y592A-K597A and A591- 
597, Rapl 

[H7-827], Rapl [358-827] and Rapl [675-827] were 
cloned into pETM-13 expression vector and expressed in 
Escherichia coli BL21(DE3) Star strain grown overnight at 
20°C after IPTG induction at [IPTG] = 0.5 mM. The 
proteins were purified using affinity chromatography 
with HisTrap column (GE Healthcare), cation exchange 
chromatography with ResourceQ columns (GE 
Healthcare), affinity chromatography HiTrap Heparine 
column (GE Healthcare) and size exclusion chromatog- 
raphy, using Superdex200 and Superdex75 column (GE 
Healthcare). Double-strand DNA were prepared using 
single strand oligonucleotides ordered at Eurogentec, 
and further purified using anion exchange chromatog- 
raphy with MonoQ 10/100 GL column (Amersham). 
Purified single-strand oligonucleotides were hybridized 
by mixing equimolar ratio of complementary DNA 
strands, heated to 96° C for 5min, and slowly cooled over- 
night. Complexes between DNA and Rapl[i_ 82 7], 
Rapl[ii7_ 82 7] or Rapl[ 358 _ 82 7] were formed by mixing 
equimolar ratio of protein/DNA, and purified by size ex- 
clusion chromatography using Superdex200 column. The 
different oligonucleotides and protein constructs used in 
the study are summarized in Supplementary Table SI. 

SAXS data collection, analysis and models construction 

Rapl [1-827], Rap l[i i7_827] 5 Rapl [358-827] alone or in 
complex with DNA were prepared in the same storage 
buffer and same concentration than for AUC experiment 
(Supplementary Methods). All data were collected on 
SWING Beamline using either the sample changer or the 
online HPLC system (43) in the case of Rapl[i_ 82 7] or 
Rapl[i_ 82 7]/DNA complex (Supplementary Figure S1A 
and B). The data were analyzed using FOXTROT 
(from SWING beamline) and PRIMUS (Primary 
Analysis and Manipulations with Small Angle Scattering 
DATA) from ATSAS 2.1 (44), from which Guinier 
and normalized Kratky plots (45) were generated 
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Chi = 0.85 Chi = 1.1 



Figure 1. SAXS analysis of Rapl and Rapl-DNA complex. (A) Schematic boxing of Rapl functional domains. (B-E) Cross-validated models from 
SAXS data, calculated with BUNCH (48), with fit and Chi residual of Rapl [358 _ 827] (B), Rapl [358 _ 827] /DNA (C), Rapl (D) and Rapl/DNA (E). The 
color code of domains is the same than in (A). Chi = l/(N- 1) x V(£((Iexp - cxlth)/sigma(lexp)) 2 ). 
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(Supplementary Figure S1A-D). From the corrected 
scattering curves, the pair-distribution functions were 
computed using GNOM (46), leading to the radius of gyr- 
ation R g , and the maximal distance D mSLX (Supplementary 
Figure SID, Table 1). The molecular weights of the par- 
ticles were derived from the extrapolated intensity at the 
origin I 0 (Table 1). In the cases of the proteins alones, 
the program DAMMIN (47) was used to generate the 
low-resolution ab initio shape envelops, which were 
superimposed and averaged using DAMAVER (44) 
(Figure S1E). To calculate the models, we used PDB 
coordinates from the RCT X-ray structure (entry 3CZ6, 
38) and from our crystal structure of DBD/DNA for the 
three constructs, and from the BRCT NMR structure 
(entry 2L42, 41) for Rapl full length and Rapl[ii 7 _ 82 7]. 
The models were calculated in an iterative manner with 
the program BUNCH (48) and cross-validated with AUC 
sedimentation coefficients (see Supplementary Methods 
for detailed procedure and Table 1). 



and K597A (GCT instead of AAG) for the double 
mutant, a 21 -bp deletion that lead to the absence of 
residues 591-597 for the deletion mutant and R580A 
(GCT instead of AGG) for the single mutant. These 
plasmids, together with a plasmid encoding wild-type 
Rapl, were introduced in yeast strain Lev559 (W303-la 
background MATa ade2-l trpl-1 ura3-l leu2-3,112 
his-3-11,15 canl-100 RAD5 rapl- (A )::KANr 
pACEl-UBRl pACEl-ROXl) (52). 

Western blot 

Exponentially growing cells in synthetic medium lacking 
tryptophan and supplemented or not with 400 uM CuS0 4 
during 6h were subjected to protein extraction. Proteins 
were extracted at 95°C in urea 8 M, SDS 2% and Tris- 
HC1 pH 7.5 100 mM followed by vortexing at 2500 r.p.m. 
with glass beads. About 3 ug of total extract were loaded 



Crystallization and structure determination 

Initial crystallization conditions were obtained at PF6 
(Institut Pasteur, Paris, France) and further optimized in 
our laboratory. The most suitable crystals were obtained 
using seating drop vapor diffusion at room temperature, 
with a concentration of Rapl[ 358 _ 8 27]/DNA complex 
between 3.3 and 5.5mg/ml, in a solution containing 
100 mM Tris-HCl pH 8.0, 20% (v/v) PEGmme-550, 
100 mM CaCl 2 and 5% (v/v) glycerol. Native diffraction 
data were collected on PX1 beamline at SOLEIL synchro- 
tron, and reduced with XDS (49). The structure was 
solved by molecular replacement with PHASER (50), 
using 1IGN pdb entry as model probe (12), and refined 
with BUSTER5 (51). Electron density and SDS page gel 
analysis revealed that Rapl[3 58 _ 8 27] construct was partly 
degraded and that in addition to the 31 -bp oligonucleo- 
tide, only the DBD moiety and part of DBD-RCT linker 
was present in the crystals. Statistics of the data collection, 
refinement and model validation are shown in Table 2. 

Directed-mutagenesis for the in vivo experiments 

RAP1 sequence (from 676 bp before the ATG to 427 bp 
after the STOP codon) was cloned into pRS414 plasmid. 
Introduced mutations were Y592A (GCT instead of TAT) 



Table 2. Statistics of the diffraction data collection, refinement and 
model validation 

Data collection 



Wavelength (A) 
Space-group 

Diffraction limits (last shell) 
Unit cells (axbxcxax(3xy) 

Emerge 

Number of unique reflections 

I/O 

Completeness 
Molecular replacement 
Refinement 
Resolution 

-^work 
^free 

Figure of merit 
Number of residues 
Number of bases 
Number of water molecules 
RMSD bond 
RMSD angles 

Average B-factor (protein, o A 2 ) 
Average B-factor (DNA, A 2 ) 
Average B-factor (residues 675-601, A 2 ) 
PDB entry 



0.98 
P2 1 2 1d 2 1 

2.95A (3.13A-2.95A) 
40.6 x 102.9 x 116.8 
90 x 90 x 90 
0.138 (0.688) 
9831 

14.32 (3.39) 

0.901 (0.688) 

LLG = 1305.76 

Buster5 

2.95 A 

0.1923 

0.2650 

0.869 

223 

62 

6 

0.0010 
1.48 
62 
100 

69 

3UKG 



Table 1. Analytical ultracentrifugation, SAXS analysis and cross-validation of SAXS models 





Construct 


R a pl [358-827] 


Rapl[ii 7 _827] 


Rapl 


Rapl [358 _ 827] /DNA 


Rapl[ii 7 -827]/DNA 


Rapl /DNA 




Th. MW (kDa) 


55 


82 


94 


67 


94 


105 




Partial specific 


0.717 


0.712 


0.710 


0.689 


0.693 


0.693 




volume (ml/g) 














AUC 


[Sample] (mg/ml) 


0.55, 1.1, 2.2 


0.82, 1.64, 3.11 


0.56, 0.94, 1.4, 
1.88, 2.81, 3.74 


0.67, 1.35, 2.02 


0.56, 0.94, 1.88, 2.81 


0.63, 1.06, 2.12, 
3.17, 4.23 




MW (kDa) 


55 


82 


94 


67 


94 


106 




Sed coef (So w 20) 


3.6 


3.8 


3.9 


4.2 


4.3 


4.2 




f/fo 


1.5 


1.9 


2.0 


1.6 


2.0 


2.0 


SAXS 


[Sample] (mg/ml) 


1, 3.2 


1.8, 2.5 


l.l a 


9 


0.7, 1.4, 2.8, 4.1 


1.2 a 




*g (A). 


32.7 


58.1 


67.4 


40.6 


71.6 


72.4 




£>max (A) 


120 


224 


260 


162 


275 


260 




Calc. sed. Coef. 


3.6 


3.8 


4.0 


4.2 


4.3 


4.2 
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on 4-12% NuPage gel (Invitrogen 1M ), and then blotted 
on a Hybond-P membrane (Amersham). Rapl was 
detected with RAP1-Y-300 rabbit polyclonal IgG from 
Santa Cruz Biotechnology. 

RESULTS 

Orientation of Rapl RCT is modulated upon DNA 
binding 

The overall structural organization of free Rapl was 
compared with that of Rapl bound to DNA by using 
analytical ultracentrifugation and small angle X-ray scat- 
tering techniques. Three different constructs were tested: 
Rapl [358-827], Rapl[ii 7 _827] and Rapl full-length (Rapl). 
Rapl [358-827] includes the DBD, the toxicity region, the 
transcription activation region and the RCT domain 
(Figure 1A), and corresponds to the required and suffi- 
cient region of Rapl for telomere associated functions 
(53). The long N-terminal moiety Rapl[i_ 35 7] is predicted 
as flexible, and was built in two steps using the construct 
Rapl [H7-827] (including BRCT, DBD and RCT domains), 
and Rapl full-length (Rapl). The molecular weights 
derived from Guinier analyzes of the SAXS data 
and the sedimentation coefficient distribution profiles 
from analytical ultracentrifugation experiments reveal no 
sample aggregation, and confirmed that these three con- 
structs of Rapl are monomeric in solution and that Rapl/ 
DNA complexes have a one to one stoechiometry 
(Supplementary Figures SIC, S2 and Table 1). 

The use of Kratky plots to detect macromolecular flexi- 
bility from SAXS data is well established and extensively 
used (54-56). The bell shapes of Kratky plots calculated 
for Rapl [358-827] free and in complex with DNA reflect 
their compact organization (Supplementary Figure SID). 
The pair distribution profile corresponds to a globular 
molecule in the case of free Rapl[ 358 -827] 5 and to an 
elongated molecule in the case of its complex with DNA 
(Supplementary Figure S1E). Average structures were 
calculated with BUNCH (48), and structures consistent 
with SAXS and AUC data were further selected 
(Supplementary Methods, Table 1). These structures 
of Rapl [358-827] and Rapl[358_827]/DNA are characterized 
by residual chi values of 0.87 and of 0.64, respectively 
(Figure IB and C). Interestingly, the RCT orientation 
with regard to the DBD is different when the molecule 
is alone or in complex with DNA; this implies a conform- 
ational adjustment of Rapl[3 58 -827] upon DNA binding. 

Unlike in the case of Rapl[ 358 -827] 5 the Kratky plots of 
both Rapl [H7-827] and Rapl, free or in complex with 
DNA, correspond to partially unstructured molecules 
(Supplementary Figure SID), and pair distribution 
profiles are characteristic of elongated molecules 
(Supplementary Figure S1E). Calculated ab initio enve- 
lopes reveal a compact shape for Rapl[ 358 -827] 5 and a pro- 
gressive elongation for Rapl[n 7 _827] and Rapl 
(Supplementary Figure S1F). Superposition of the three 
envelops highlights structural conservation of a compact 
region that could be attributed to region [358-827] 
(Supplementary Figure S1F). The frictional ratios 
measured by AUC arise from 1.5 to 1.9 and 2.0 in 



Rapl [358-827], Rapl[ii 7 _827] and Rapl, respectively, 
showing the transition from a compact structure for 
Rapl [358-827] to elongated structures for Rapl[ii 7 _ 8 27] and 
Rapl (Table 1, Supplementary Figure S2). The average 
structures of Rapl[ii 7 _ 8 27] and of Rapl, free or in 
complex with DNA, remain compact in region 358-827, 
the N-terminal half is extended in all models 
(Supplementary Figures S3 A and B, S1D-E), and the 
overall shapes of the molecules free and in complex 
with DNA are similar (Figure S3C-D). We calculated the- 
oretical SAXS curves using region [358-827] from 
Rapl [H7-827] and Rapl models free or in complex with 
DNA, and compared them respectively with experimental 
Rapl [358-827] or Rapl[358_827]/DNA SAXS curve. We 
observed a better fit between theoretical and experimental 
curves in the case of DNA complexes, although super- 
position of region [358-827] from free Rapl[3 58 -827] ? 
Rapl[H7_827] and Rapl reveals similar overall shapes 
(Supplementary Figure S4A and B). Superposition of 
region 358-827 on the DBD domain from the three 
Rapl constructs shows that the relative orientation of 
RCT with regard to DBD fluctuates when the protein is 
free, whereas it is confined to a limited area in the DNA 
complex (Supplementary Figure S4C-D). Moreover, the 
orientation of the RCT domain does not depend on the 
presence of the N-terminal moiety of the molecule. 

In summary, comparison of the average structures from 
our three constructs highlights two important points: (i) 
the compact conformation of region 358-827 does not 
depend on the presence of the N-terminal moiety; (ii) the 
interaction with DNA is associated with a reorientation of 
the RCT and a decreased flexibility. 

Structural determinants of Rapl organization upon 
DNA binding 

The reorientation of the RCT domain with regard to the 
DBD upon DNA binding suggests that the overall con- 
formation of Rapl is driven by important structural deter- 
minants located in the DBD-DNA moiety. The available 
high-resolution crystal structure of DBD in complex with 
a 19 bp double-stranded DNA shows that the two Myb- 
domains are tandemly arranged on the DNA so that each 
is aligned to and makes similar contacts with the GGTGT 
sequences of Rapl binding site (12,57,58). However, due 
to crystal packing interactions, the available structure 
leaves uncertainty concerning the conformation of the 
C-terminal loop of DBD (12), which limits our knowledge 
about the RCT orientation with respect to the DBD. 
To clarify this information, we solved the X-ray structure 
of Rapl dbd m complex with a 31 bp double stranded 
DNA oligonucleotide (Table 2). The oligonucleotide 
length enlarges the binding interface, and reduces the 
crystal packing constraints that tighten the DNA con- 
formation. The superimposition of our structure with 
the available Rapl D BD/DNA complex coordinates (PDB 
entry 1IGN, 12) leads to a rms deviation of 1.77 A. As 
expected, Rapl dbd exclusively binds the two hemi-sites 
spaced by three nucleotides; the additional DNA region 
remains free of protein (Figure 2 A and B). 
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Figure 2. Structural determinants of Rapl-DNA interaction. (A) Overall structure of Rapl-DBD in complex with DNA, with lsigma level electron 
density map in blue around the region corresponding to the wrapping loop (residues 565-601). H1-H7 labels highlights alpha-helices in sequential 
order. The blue and red arrows indicate the respective directions of N- and C-terminal ends. (B) Zoom on Rapl wrapping loop from Lys575 to 
Asn601 (in green). 



The main novelty of our structure lies in the quality of 
the definition of loop [565-601]. The 1IGN structure 
lacks three protein regions: loops [482-512], [565-571] 
and [579-586], and ends at position 594. In our structure, 
the loop corresponding to [482-512] is still disordered 
between residues 480-498. The [565-601] region (seven 
residues longer than in 1IGN) is well defined, and wraps 
around the DNA molecule. It runs along DNA major 
groove with no secondary structure element, and brings 
the N and C-terminus extremities nearby, although 
they are orientated in opposite directions (Figure 2A). 
In addition, this wrapping loop is deeply engaged in 
protein/protein and protein/DNA interactions, with a 
total buried surface of 1508 A 2 , corresponding to 37% 
of its total surface. The interaction with DNA buries 
2374 A 2 of the DBD, among which 642 A comes from 
the wrapping loop. The wrapping loop starts at region 
565-571, lengthens helix 7 by one helical turn, and 
forms a tight turn at residues 565-569 that directs the 
loop back toward DNA major groove (Supplementary 
Figure S5A). The interaction between the wrapping loop 
and DNA involves segments 575-583 and 589-597. 
Residues 575-583 interact with Adel9, Cyt20, Ade21 



from C-rich strand and Gua8 from G-rich strand, and 
Arg580 is deeply plugged in the major groove 
(Supplementary Figure S5B). Region 589-597 interacts 
with Cytl6 and Cytl7 from C-rich strand, and residues 
Tyr592 and Lys 597 interact with helix 2 from Mybl. 
These Mybl /wrapping-loop interactions close the clamp 
formed by the wrapping loop (Figure 3A). Finally, the 
N-terminal and C-terminal extremities of DBD are 
oriented in opposite direction, in agreement with the ob- 
servation provided by our SAXS analysis. 

Rapl wrapping loop is important for its functional 
integrity and its binding to DNA 

Our SAXS analysis together with the crystal structure 
implies that the wrapping of Rapl around DNA has im- 
portant consequences regarding its overall conformation, 
and its adjustment upon DNA binding. The physiological 
relevance of the wrapping loop was addressed with three 
Rapl alleles mutated within the loop: Arg580 from the 
first segment was modified into alanine (ra/?/-R580A), 
residues Tyr592 and Lys597 that lock the clamp were 
replaced by alanines (rapl-Y592A; K597A), and the 
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Figure 3. Functional importance of Rap 1 wrapping loop. (A) Cartoon representation of C-terminal clamp region colored from conservation analysis 
scores. Residues involved in C-terminal clamp are represented in sticks. Red and blue arrows highlight residues of the wrapping loop (red) that 
interact with residues of Mybl (blue). (B) Western blot comparing Rapl amount in the different strains. Rapl level in rapl-(A ) cells complemented 
with RAP1 (WT), rapl-Y592A;K597A, rapl-A591-597 and rapl-R580A alleles. The asterisk indicates the position of the tagged Rapl degron, which 
disappears upon Cu 2+ addition. (C) Ability of the rapl alleles to complement Rapl loss. The 10-fold serial dilution of exponentially growing cells was 
spotted on synthetic medium plates with or without 400 mM CuS0 4 . The plates were photographed after 3 days of growth at 30° C. (D) ITC titration 
of Rapl, Rapl [Y5 92A-K597A] or Rapl 

A[59i-597] by DNA. (E) EMSA experiments using 8 |iM DNA in the presence or not of 4|iM of Rapl, 
Rapl[Y592A-K597A] or Rap 1 A[59i-597]- Cplx indicates the position of DNA/protein complex; DNA indicates free DNA. 



whole region 591-597 was deleted (rapl-A59 1-597). These 
constructs were introduced in Lev559, a strain encoding a 
Rapl degron protein that will be degraded upon copper 
addition (14), allowing in vivo analysis of mutated Rapl 
proteins. As shown in Figure 3B, rapl-R580A protein 
level is identical to wild-type level but rapl-Y592A- 
K597A and rapl-A591-597 levels are slightly increased. 
All strains display normal growth when rapl degron is not 
induced. Upon Cu 2+ addition, uncomplemented cells die 
whereas cells complemented with wild-type RAP1 or with 
the rapl-R580A allele grow normally (Figure 3C). 
Interestingly, the rapl-Y592A-K597A and rapl-A591- 
597 alleles fail to fully compensate endogenous Rapl 
loss, the rapl- A59 1-597 mutant displaying a more severe 
growth defect. In the same way, ITC titration of Rapl 
wild type or mutants Rapl[Y592A-K597A] and 
RaplA[591-597] by DNA reveals a lower affinity by a 
factor of 2 or 5, respectively for each mutant as 
compared to the wild type (Figure 3D, Supplementary 
Table S3). Native gel migration profiles are similar for 



both wild type and mutated Rapl /DNA complexes, 
although the size and flexibility of the molecule suggest 
that changes in domain orientation remains possible 
(Figure 3E). These analyses indicate that the region 
[591-597] as well as residues 592 and/or 597 are important 
if not essential for Rapl function, and highlight the im- 
portance of DBD wrapping loop and its clamp through 
interactions with Mybl. 

Orientation of C-terminal domain within Rapl /DNA 
complex is compatible with its interaction with functional 
partners Rif2 and Sir3 

Two significant features come from our analysis: (i) the re- 
orientation of the RCT domain with regard to the 
DBD upon DNA binding, (ii) the location of Rapl N- 
and C-terminal parts on opposite sides of the DNA 
(Figure IE). Numerous studies have shown that RCT 
domain plays a crucial role in Rapl partners recruitment 
(6,15,29-35). We therefore checked if the conformation of 
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Rapl/DNA complex and the orientation of RCT was 
compatible with its interaction with functional partners 
in solution. We performed NMR titrations of 
15 N-labeled Rapl[ 67 5_ 8 27] with a Sir3[ 450 -487] interacting 
peptide designed on the basis of previous double-hybrid 
experiments (59) and Rif2 full length (Supplementary 
Methods, Supplementary Figure S6). After assigning the 
NMR spectra of Rapl[ 675 _ 82 7] (Supplementary Methods, 
Figure S6A), we identified residues whose peaks were 
affected upon titration with Sir3[ 450 -487] an d Rif2. The re- 
sulting surfaces are consistent with previous studies 
(38,40), and are partly overlapping (Supplementary 
Figure S7A). In addition, these Sir 3 and Rif2 interacting 
surfaces on Rapl RCT are still accessible when positioned 
onto the average structure of Rapl full length in complex 
with DNA (Supplementary Figure S7B). 

DISCUSSION 

Interaction of Rapl with DNA is associated to structural 
adjustments 

The integrated structural biology study presented here 
allows for the first time to build a complete molecular 
image of the average architecture of Rapl in solution 
free and in complex with DNA, despite its modular or- 
ganization and flexible character. In the absence of DNA 
the molecule is elongated and partially disordered. The 
N-terminal part is particularly extended and poorly 
organized, whereas the DBD-RCT is more compact. 
The combined analysis of the three different constructs 
confirms the modular character of the molecule and also 
shows that the relative orientation of the globular 
domains fluctuates one with respect to the others. 
Comparison of our results with available electron micros- 
copy images of Rapl is of particular interest since our 
approach describes the average extended conformation, 
whereas high-resolution microscopy highlights the more 
frequent pseudo-ring (13,39) or V-shape (22) conform- 
ations. NMR titrations, ITC and gel filtration experiments 
confirmed that there is no preferential interaction between 
the N-terminal part of Rapl and the RCT domain or 
DBD-RCT region that could stabilize a pseudo-ring 
structure of Rapl (details in Supplementary Methods 
and Supplementary Figure S8). The extended and 
modular arrangement of free Rapl, as well as the accessi- 
bility of the Myb domains, is favorable to its interaction 
with DNA. However, modeling a DNA/Myb interaction 
using the domain organization of free Rapl leads to 
sterical hindrance between DNA and RCT domain 
(Figure 4A). This shows that the average domains organ- 
ization of Rapl needs to adjust upon DNA binding, in 
agreement with the conformation observed for Rapl in 
complex with DNA (Figures IE and Supplementary 
Figure S3B). The conformational adjustment of Rapl 
upon DNA binding leads to an orientation of the N- 
and C-terminal parts of the molecule in opposite direc- 
tions almost perpendicular to the DNA axis (Figure IE). 
Our crystal structure of DBD/DNA reveals that such con- 
formational adjustment is probably driven by a loop 
located at position 575-601 of DBD that wraps DNA 



molecule (Figure 2A). We observed a reorientation of 
the RCT relatively to the DBD due to DNA binding, 
and the orientation of the RCT is more constrained with 
regard to the DBD. Furthermore, NMR titration of 
15 N-labeled Rapl[ 67 5_ 8 27] with the full length protein 
Rif2, and a Sir3[ 450 -483] peptide leads to interacting 
surfaces, which remain solvent-exposed in our average 
structure of Rapl/DNA complex and thus accessible to 
Rif2 and Sir3 (Figures S6 and S7). Therefore, we 
propose that the orientation of the RCT in Rapl/DNA 
complex may facilitate its interaction with its functional 
partners, and the corresponding telomeric functions. 

Remodeling of Rapl involves its wrapping around DNA 

The wrapping loop located at position 575-601 of DBD is 
involved in extensive interactions with DNA, and in- 
creases the surface involved in DNA binding by 25%. In 
addition, its interaction with the first Myb of the DBD 
clamps the molecule around DNA through two 
residues (Tyr592 and Lys597) located at the tip of the 
loop (Figures 2 A and 3A). In vivo experiments using 
yeast strains expressing mutants of Tyr592 and Lys597 
or deletion of 591-597 segment highlight the functional 
significance of these residues (Figure 3B and C). 
Importance of region 583-596 in Rapl DNA binding 
was previously observed using band shift assays (60). In 
this publication, Kingsman's group observed DNA 
binding with Rapl region 302-596, but not with region 
302-583, and concluded that the DNA binding domain 
ends within segment 583-596. In our study, we observe 
that mutation of Rapl at Tyr592 and Lys597, or 
deletion of 591-597 segment affect its affinity for DNA 
by a factor of two and five respectively, although DNA 
interaction is not fully abolished (Figure 3D and E). 
Altogether, our results suggest that (i) the wrapping, 
(ii) the clamp of the molecule at a specific position on 
DNA molecule and/or (iii) the orientation of N- and 
C-terminal part in opposite direction are essential for 
Rapl functional integrity. In fact, in the case of 
telomere protection, chromatin opening, or interaction 
with nucleosomal binding sites, the tight wrapping of 
Rapl dbd is probably sufficient since these functions can 
occur in the absence of N- and C-terminal domains of the 
molecule (9,26,29,36). 

Toward a molecular model of telomere assembly 

Combination of NMR results, crystal structure and SAXS 
analysis allows us to build the overall architecture of Rapl 
in interaction or not with double-stranded DNA. These 
results highlight the modular character of Rapl, its struc- 
tural plasticity that may facilitate its binding to DNA, and 
a complete and tight wrapping of the molecule around 
DNA that constraints the orientation of the C-terminal 
fragment and may favor its interaction with its functional 
partners (Figure 4A). 

In a telomeric context, Rapl binds independently to 
each recognition site (11, 13). In the same way, our inter- 
action studies show neither oligomerization of Rapl 
(Supplementary Figure S2, Table 1), nor interdomain 
interaction. In addition, the extended conformation of 
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Figure 4. Architecture of telomere first shell macro-assembly. (A) Complete architecture of free Rapl, model of free Rapl in interaction with DNA 
and Rapl/DNA complex in interaction with a telomeric DNA, with the RCT surface involved in functional partners interactions, based on the 
combination of SAXS analysis, crystal structure and NMR titration. DNA is represented in surface mode colored in grey and telomere hemi-site in 
purple, BRCT in orange, DBD in green and RCT in red (with overall surface involved in Rif2 and Sir3 interaction in blue). The path of the 
wrapping loop is indicated with a purple arrow. (B) Modeling of the interaction of a second Rapl molecule along telomere. (C) Reconstitution of a 
four Rapl sites telomere along or in front of DNA fiber. Bound Rapl molecules are represented in cartoon mode with different colors. Rapl sites in 
violet are spaced every 18 bp. 



the molecule bound to DNA, with its axis of elongation 
almost perpendicular to the axis of the DNA molecule, 
leaves accessible the nearby sites, and is compatible 
with the binding and remodeling of a second molecule 
(Figure 4B). On the basis of these data, we reconstituted 
a molecular model of an 80 base pairs telomere with its 
first shell bound Rapl molecules every 18 base pairs (11) 



(Figure 4C). We observe that the 18 bp spacing leads to a 
high molecular density. Because of the use of a short DNA 
in our study, the limit in detectable affinities, and the tight 
binding of Rapl-DBD at DNA, one cannot exclude that 
transient interactions might occur between adjacent mol- 
ecules. The proximity between Rapl molecules may also 
favor particular macromolecular assemblies at the second 
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shell of telomere structure that involve Rapl functional 
partners. The formation of such molecular scaffolding 
now needs to be documented at a structural level. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Methods, Supplementary Figures 
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